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Description 

[0001] This invention relates to a method of detecting and locating stator faults in dynamoelectric machines during 
operation thereof, and to an apparatus for implementing the method. Specifically, the invention is concerned with the 

s detection and location of a shorted turn of the stator winding by monitoring the axial leakage flux. 

[0002] Studies on reliability and failure mechanisms of dynamoelectric machines with multiturn stator coils indicate 
that stator insulation breakdown is one of the major causes of failure. Additionally there is evidence that the majority 
of stator winding failures in such machines result from turn insulation breakdown. In some cases where there has been 
minimal or no core damage due to the fault, it is possible to isolate the damaged coil from the rest of the winding and 

10 restart the machine. Such an operation may necessitate isolation of other coils in order to retain balanced operation 
of the machine. In this way the machine can be returned to service quickly, albeit with degraded performance until a 
more complete repair or winding replacement can be scheduled. Unfortunately location of the failed coil is difficult and 
time consuming. 

[0003] IEEE Proc. B. El. Power Appl., Vol. 133, no. 3, pages 142 to 148 discloses a method of monitoring electric 
75 drives in order to pre-empt occurrences of failures and to give an early indication of the onset of possible failures. Fault 
signatures are used to identify the type of fault which may occur and to give an early indication of the onset of a possible 
malfunction. Specifically, a possible malfunction is detected by measuring the harmonic components in axial leakage 
flux and comparing the measured harmonics with known fault signatures. 

[0004] The object of the present invention is to provide a method of, and a device for, detecting the occurrence of a 
20 shorted turn and indicating its location in the stator winding of the machine while the machine is still running. 

[0005] According to a first aspect of the present invention, there is provided a method of detecting and locating stator 
faults in a dynamoelectric machine having a rotor and a stator with multiturn stator coils by monitoring changes in axial 
leakage flux during operation of the machine, comprising: 

25 identifying the frequencies of selected harmonic components of axial leakage flux that are to be significantly af- 

fected by the occurrence of a stator fault, 

detecting the axial leakage flux at each of at least four positions, including a datum position, distributed symmet- 
rically around the axis of the rotor adjacent to an end winding of the stator, 

deriving from the detected flux at each of said positions a first signal having a value corresponding to the magnitude 
30 of the detected flux, 

storing the values of said first signals derived prior to the occurrence of a stator fault, 

deriving from the detected axial leakage flux a further signal having a value corresponding to the magnitude of 
said selected harmonic components, 

monitoring said further signal and detecting a change therein thereby detecting the occurrence of a stator fault, 
35 storing the values of said first signals derived after the occurrence of the stator fault, and 

comparing the values of said first signals stored before and after the occurrence of a stator fault and storing the 
differences therebetween, 

characterized in that the method further comprises: 

40 

determining from said comparison a value representing the angular position of the stator fault with respect to the 
datum position by: 

(i) storing a value R-, representing the radial distance of said four positions from the axis of the stator, 
45 (ii) storing a value R 2 representing the mean radius of the stator end winding, 

(iii) computing values B-, and B 2 corresponding respectively to the moduli of said differences between the 
stored values of said first signals derived from the detected flux at one diametrically opposed pair of said flux 
detecting positions, 

(iv) computing values B 3 and B 4 corresponding respectively to the moduli of said differences between the 
50 stored values of said first signals derived from the detected flux at the other diametrically opposed pair of said 

flux detecting positions, 

(v) computing values and \j/ 2 , where 

55 ± ¥1 = cos'" l [0.5(R 1 /R 2 +R 2 /R 1 ).(B r B 2 )/(B 1 -»-B 2 )] 



± V 2 =cos" 1 [0.5(R 1 /R 2+ R 2 /R 1 ).(B 3 -B 4 )/(B 3+ B 4 )], 
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(vi) determining therefrom the angular position of the stator fault \j/ 0 with respect to said datum position from 
the relationship 

5 Vo = Vi = V 2 + 90 ° 

and (vli) displaying the value \j/ 0 . 

[0006] The occurrence of a stator fault representing a shorted turn is thus detected by monitoring selected axial flux 
10 harmonics which are sensitive to the occurrence. The angular position of the fault is computed from stored values, 
namely machine parameters which are initially known and measured values of the flux after detection of the fault. 
[0007] According to a second aspect of the present invention, there is provided an apparatus for detecting and lo- 
cating stator faults in an induction motor by monitoring changes in axial leakage flux during operation of the motor, 
comprising: 

15 

means for detecting the axial leakage flux during operation of the motor, said means including means for deriving 
a plurality of first signals each having a value corresponding to the detected flux at a respective one of at least 
four positions, including a datum position, distributed around the axis of the motor, 

means for deriving from the detected flux a further signal corresponding to the magnitude of the detected flux at 
20 selected harmonic frequencies, 

first data processing means for monitoring said further signal and for detecting the occurrence of a stator fault by 
a change in the value thereof, and 

second data processing means responsive to said first data processing means, said second data processing means 
including: 

25 

(a) means for storing said first signals prior to the occurrence of a stator fault, and 

(b) means responsive to the detection of a stator fault for storing said first signals after the occurrence of the 
fault, 

30 characterized in that said second data processing means further includes: 

(c) means for computing from said stored values of said first signals the angular position of the stator fault 
with respect to said datum position, 

(d) said means for computing including: 

35 

(d1) means for storing a value representing the radial distance of said four positions from the axis of 
the stator, 

(62) means for storing a value R 2 representing the mean radius of the stator end winding, 
(d3) means for computing values B 1 and B 2 corresponding respectively to the moduli of said differences 
40 between the stored values of said first signals derived from the detected flux at one diametrically opposed 

pair of said flux detecting positions, 

(d4) means for computing values B3 and B 4 corresponding respectively to the moduli of said differences 
between the stored values of said first signals derived from the detected flux at the other diametrically 
opposed pair of said flux detecting positions, 
45 (d5) means for computing values and \j/ 2 , where 

± \|f 1 = cos" 1 [0.5(R 1 /R 2 +R 2 /R 1 )»(B 1 -B 2) /(B 1 +B 2 )] 

so - 

± y 2 = cos' 1 [0.5(R l /R 2+ R 2 /R 1 )-(B 3 -B 4 )/(B 3+ B 4 )], 

(d6) means for determining therefrom the angular position of the stator fault \\f 0 with respect to said datum 
position from the relationship 

55 
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and (d7) displaying the value y 0 . 

[0008] The preferred method according to the invention, and an apparatus therefor, will now be described by way of 
example with reference to the accompanying drawings. In the drawings: 

5 

Figure 1 is a graph showing the waveform of detected magnetomotive force (MMF) due to a shorted stator turn; 
Figures 2 and 3 are diagrams illustrating the geometric considerations required in magnetic field calculation; 
10 Figure 4 is a diagram illustrating the geometry of sensors used to detect axial leakage flux; 

Figure 5 is a partially diagrammatic illustration of an apparatus used to monitor axial leakage flux; 
Figure 6 is a diagram illustrating the general algorithm used to detect and locate a stator fault; and 

15 

Figure 7 is a diagram illustrating the detailed algorithm used to locate the stator fault. 
Basic Principles 

20 [0009] In the case of an ideal electrical machine there is zero net axial magnetic leakage flux because the stator and 
rotor currents should exactly cancel one another However, in practice, due to nonuniformities in materials and con- 
struction methods, such a situation does not occur. Hence associated with all rotating machines there is a small, but 
measurable, axial leakage flux. This magnetic flux can be measured by the simple expedient of placing a search coil, 
or other appropriate magnetic flux sensing device, in close proximity to the end region of the machine. The output from 

25 such a device, when analyzed in the frequency domain, is manifest as a series of harmonics of varying magnitude. 
From basic electromagnetic considerations of electrical machines it is possible to derive expressions which permit the 
harmonic components inherent in the axial leakage flux to be predicted. The equations, which are simple algebraic 
expressions, only require knowledge of the number of poles and rotational speed of the machine, in order to calculate 
all possible harmonics for a particular machine. 

30 [0010] The occurrence of a stator fault in an induction motor, for example, will inevitably result in a change in the air 
gap space harmonic distribution. These space harmonics cannot be detected directly by a search coil. However, the 
search coil can detect the time harmonics of the axial flux. Thus it is necessary to derive the relationship between the 
space and time harmonics in order to correctly interpret the frequency spectrum obtained from the search coil. For the 
purpose of this work, considering only the supply fundamental and the third harmonic component due to saturation, 

35 the space harmonic distribution of MMF due to a balanced, full pitched, three phase winding fed from a balanced supply 
frequency, oe>, is given by [1], 

m = 0.955N 1 [k w1 cos(cot-p0)+O.2k w5 cos(Q)t+5pe) 

40 

-0. 1 4k w7 cos(cot-7pe)+0.09k w1 1 cos(cot+1 1 p8)-. . .] (1 ) 

where 

45 k^ is the n th winding factor 

p is the number of pole pairs 
9 is the angular displacement from the stator datum 

[0011] This represents a rotating set of harmonics of order 6n±1 which can be simplified to the corresponding air 
50 gap fluxes, 

B Q = B 1 cos((ot-pe)+B 5 cos(cot+5p0)-B 7 cos(cot-7pe) 
55 +8^006(0)1+1 1p0) (2) 



where B n are the spatial harmonic fluxes. 
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[0012] This expression is in the stator frame of reference. Consequently, because the shaft flux of the rotor is of 
interest, it is necessary to refer equation (2) to the rotor reference frame. Consider the situation in which p is the angular 
displacement between the rotor and stator datum positions, and a is defined to be the angular displacement from the 
rotor datum. Then 9 = a + p. 
5 [001 3] If the angular rotor speed is co r , then, 

0 = a+o> r t (3) 

10 now using the normal expression for the slip of the motor, i.e. s = (Gva> r )/oo s , where gd s , the synchronous speed, 

= oVp, 

co r = G)(1-s)/p (4) 

15 

[0014] Now, the general term of equation (2) is, 

B ns = B n C0S ( G)t±n P e ) < 5 ) 

20 

substituting equations (3) and (4) into (5) leads to, 

B ns = B n cos[(1±n(1-s))(ot±npcc] (6) 

25 

[0015] Expanding the expression for the first few terms gives, 



30 



B s = B 1 cos(scot-pa)+B 5 cos[(6-5s)(ot+5pa] 
-B 7 cos[7s-6)cot-7pa]+B 1 1 cos[(1 2-1 1 s) 
+ 11 pa]- 



(7) 



[0016] Equation (7) gives the frequency components of the currents that are induced in the rotor due to the air gap 
space harmonics of a balanced winding and supply. In addition to these harmonics, the fundamental of the supply 
frequency will also appear in the axial flux spectrum. The presence of additional, higher order, harmonics can be 
accounted for by using the term nco rather than co. 

40 7] Having derived the harmonic components of axial flux due to asymmetries in a "healthy" machine, the situation 

in which an interturn short circuit is introduced is now considered. The effect of the interturn fault is to remove a turn 
from the stator winding. This will have a small, but finite, effect on the main air gap flux distribution. In addition, an emf 
will be induced in the shorted turn which will result in a current flow limited only by the self impedance of the fault. This 
impedance essentially determines the transition time between turn and groundwall insulation failure. 

45 [0018] The fault current due to the shorted turn is the source of an additional MMF pulse. This MMF pulse also has 
a space harmonic distribution which is superposed on the main field distribution. From previous considerations, this 
will lead to a change in the time harmonics observed in the leakage field. The changes to be expected can also be 
predicted mathematically, as follows: 

[0019] Simple consideration of the MMF distribution due to an interturn short circuit leads to the characteristic illus- 
SQ trated in Figure 1. This is the case for a four pole machine. The analysis to be outlined below is for the general case 
of a 2p pole machine. 

[0020] Fourier analysis of the MMF waveform, illustrated in Figure 1 , shows that it contains all harmonics except the 
fourths, i.e., 

S5 B s = 0.5ZB n cos(cot±n9),n * 4m, all m (8) 



[0021] For the general case, the corresponding waveform would have a mark-space ratio of 1 :(2p-1 ) causing every 
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2p th harmonic to be absent. Hence the time harmonics produced by the rotor are given by, 

B s = 0.5ZB n cos[(1 ±n(1 -s)/p)cot±na], n * 2pm, (9) 

[0022] Adding in supply time harmonics of higher order k, leads to the completely general expression, 

B s = 0.522;B n cos[(k±n(1 -s)/p)cot±na], n * 2pm, (1 0) 

[0023] Although this is a large series, only the lower order harmonics are significant. The key element of this expres- 
sion is, 

75 [k±n(1-s)/p]f 1 

for k= 1,3 and n = 1,2,3, ,(2p-1) (11) 

[0024] The net term in the argument of (10) causes the components defined above to beat at the slip frequency of 
20 the rotor current. 

[0025] Thus the following procedure can be employed to identify the stator fault: 

1 . Calculate air gap space harmonics due to balanced supply and winding conditions. 

2. Derive associated time harmonic currents (fluxes) in rotor and stator reference frames. 



25 



30 



3. Refer all time harmonic components to the stator frame, since the search coil is stationary. These are the com- 
ponents to be expected in the "healthy 0 machine. 

4. Calculate additional space harmonics injected by the occurrence of an interturn short circuit. 

5. Relate these to additional time harmonics in the spectrum of the axial field. 

6. Look for changes in these components, which can be predicted for a given machine, to indicate the presence 
35 of an interturn fault. 

[0026] Having determined that a stator fault condition exists, the next task is to locate the position of the fault within 
the stator winding. The introduction of a shorted turn in the stator winding produces an asymmetry in the magnetic field 
in the endwinding of the machine. This is because the short circuit current flowing in the fault is not balanced by an 
equal current flowing in the corresponding phase belt which is situated diametrically opposite the faulted turn. Thus, 
by measuring the asymmetry in the endwinding magnetic field the position of the faulted coil within the stator winding 
can be determined. 

[0027] The technique relies upon the use of an array of at least four search coils, or other appropriate sensors, 
distributed axisymmetrically about the drive shaft of the motor. These coils provide a local measure of the magnetic 
field in the end region of the machine. By triangulating the outputs from the search coils, it is possible to determine the 
location of the stator coil containing the interturn short circuit. In order to accomplish this function it is necessary to 
develop an expression for the field at any point on the circumference of the circle passing through the centres of the 
search coils due to the current flowing in the end winding of an arbitrarily positioned turn of the stator winding. By 
invoking the Biot-Savart law and with reference to Figure 2, one can deduce 



40 



45 



50 



dB = {[l q ldjxa z )/4;ir 2 (12) 



where 

55 



B is the magnetic field strength vector 
H 0 is the permeability of free space 



5 



10 
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a z is the unit vector in the z direction 
[0028] From Figure 2, and rewriting equation (12) 

dB= (^ 0 ldlsin<t>)/47tr 2 (13) 

but since dlsin<l> = rd6, 

dB=(u Id6)/4nr (14) 



[0029] Figure 3 illustrates the general case in which the shorted turn is displaced an angular amount from an 
arbitrary datum position. In this case R 2 is the mean radius of the end winding, and R 1 is the radius of the circle on 
15 which the field is to be calculated. The circles are assumed to be coplanar and concentric. The length of arc of the coil 
is L and the field measurement point is displaced an angle 8 from the datum. Using the cosine rule, 



20 



25 



30 



35 



50 



r 2 = R 1 2 +R 2 2 -2R 1 R 2 cos(\j/+5) (1 5) 



[0030] If R 2 is very much greater than R v then L = r6, thus 



6 = L/[R 1 2 +R 2 2 -2R 1 R 2 cos(\j/+5)] 0 - 5 (16) 



[0031] From equations (14) and (16), the value of B at the field point is given by, 



B= ii 0 \U[4n{R*+R 2 2 -2n,R 2 cos(y+b))] (17) 

[0032] By fixing the position of a number of coils with respect to the datum, and arbitrarily assigning a reference coil 
at 8 = 0, the value of B at each coil can be found as a function of y. From the measured values of B, \j/ can be found 
and consequently the position of the faulted turn. 

[0033] In order to optimize the detection scheme, it is necessary to determine the minimum number of coils required 
for simple and reliable fault location. Consider the situation in which two coils are spaced 180° apart. Each coil will 
have an induced emf, proportional to the flux linking it. With reference to Figure 4, and taking the coil at position 1 to 
be the datum, from equation (17), 

40 B 1 = k/(a-bcos \j/) (18) 

and 

45 B 2 = k/(a-bcos(\|/+180°)) (19) 

thus 

B 2 = k/(a+bcos\|/) (20) 

where 



55 



k = ^ 0 IU47t 
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a = (R 1 2 + R 2 2 ) 
b = 2R 1 R 2 

[0034] From equations (18) and (20), and some rearrangement of the terms, 

cos ¥ = [a(B r B 2 )y[b(B 1+ B 2 )] (21) 

[0035] This expressions locates the fault to within +\|/ because, 

± V = cos" 1 [ 0.5(R 1 /R 2 +R 2 /R 1 ){(B 1 -B 2 )/(B 1 +B 2 )}] (22) 

[0036] In practice the measured values of the emf induced in the coils can be used because the B values are ratioed 
in equation (22). 

[0037] Consequently the values used when attempting fault location should be the modulus of the difference between 
the "healthy" and "faulted" conditions. In order to locate uniquely the fault, a second set of coils, 3 and 4, is required, 
as shown in Figure 4. For convenience they are displaced at right angles to the axis of symmetry of coils 1 and 2. 
However, the angular displacement of the coils is not critical. Repeating the location procedure in the same fashion as 
above with the second set of coils, will yield +vj/ p2 , which differentiates it from ±y p1 , found from coils 1 and 2. If the 
higher magnitude coil signal is chosen as the point from which to measure, then \j/ will always fall within the range of 
0° to 90°. In addition, only one quadrant will contain an angular position identified using both sets of coils. This is the 
approximate location of the fault. 

[0038] In order to clarify the procedure for locating the position of a faulted coil within the stator winding, reference 
should be made to Figure 4. For the purpose of the location algorithm, the coils are grouped in pairs. In this case coils 
1 and 2 form one pair, and coils 3 and 4 form the other pair. 

[0039] Initially, at a point in time at which the motor is considered to be healthy, the induced emf value on each of 
the coils is measured and recorded, in the event of an interturn short circuit occurring, the emf values on each of the 
coils are again recorded. These latter values are subtracted from the initial values and the modulus of the resultant is 
taken. Subsequently, the four positive numbers thus derived, i.e. B-,, B 2 , B 3 and B 4 , are input into the fault location 
algorithm. This operation results in two expressions, 

± ¥p1 = cos" 1 [0.5(R 1 /R 2+ R 2 /R l )((B 1 -B 2 )/(B 1+ B 2 ))] 

and, 

± ¥p2 = cos _1 [0.5(R 1 /R 2+ R 2 /R 1 )((B 3 -B 4 )/(B 3+ B 4 ))] 

[0040] The resultants from these expressions are angles. The reference points are taken as the coils, in each pair, 
upon which the largest emf value was recorded. Thus, if for example coil 1 had the largest reading between coils 1 
and 2, one would say that the position of the fault was ±\j/p-, from coil 1 . Similarly, if the reading on coil 4 was the larger 
of the readings on coils 3 and 4, the fault would be located ±\j/ p2 from coil 4. Hence, it is a simple step to draw the 
conclusion that, in this case, the the fault must lie in the quadrant between coils 1 and 4 at the angular position \j/ pl -y p2 /2. 
[0041] To summarize.the detection and location scheme is as follows: 

from machine parameters calculate the axial leakage flux components under normal conditions, 

calculate which of the harmonics will change as the result of an interturn short circuit search for these components 
in the axial leakage field of the machine, and 

upon fault indication apply the fault location technique described above. 
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Implementation of the Method 

[0042] Figure 5 illustrates a preferred apparatus used to detect and locate the occurrence of a shorted turn of the 
stator winding 1 1 of an induction motor 1 0 by monitoring changes in the axial leakage flux during operation of the motor. 
5 The apparatus comprises two subsystems, namely means for sensing the axial leakage flux and means for processing 
the signals derived therefrom. A further subsystem, which may be necessary in some cases, is some form of speed 
measurement in order to gain a measure of the slip of the motor when loaded. 

[0043] The flux sensors shown in Figure 5 are air cored coils 13, which are inexpensive and easily fabricated. The 
preferred location of the search coils is axisymmetric with the drive shaft of the motor, either inside or outside the motor 

10 casing. In the illustrated embodiment there are four search coils 1 3 distributed symmetrically around the axis of the 
motor adjacent to an end winding of the stator 11 . In principle there could be more than four search coils, but the use 
of four search coils arranged as diametrically opposed pairs simplifies the signal processing procedures and circuitry. 
[0044] Speed measurement, if required, can be accomplished in a number of ways. In some cases, speed information 
is available at the motor control centre in the plant. However, it is possible that some type of transduction system may 

15 be necessary. A measure of the rotational speed of the rotor is needed in order to calculate the slip of the motor. The 
slip is required because the harmonics of the axial leakage flux are dependent on the slip of the motor. Rotational 
speed can easily be obtained by measuring the time required for a reference point on the drive shaft to complete one 
revolution. In practice this can be done by painting a stripe on the drive shaft and using a device such as a photo 
transistor to generate an electrical pulse each time the stripe passes under the photo transistor. 

20 [0045] It is necessary to condition the voltage signals derived from the search coils, which are normally of only a few 
millivolts. In order to bring these signals to more manageable levels they are amplified by preamplifiers 1 4. Furthermore, 
in order to remove the power frequency component, which is relatively high and likely to saturate subsequent amplifier 
stages, the signals are filtered by power frequency notch filters 15. Alternatively, very wide dynamic range amplifiers 
may be used. The four analogue signals, after conditioning, are digitized to produce respective "first 0 signals having 

25 magnitudes corresponding to the harmonic flux detected by the four search coils 13. These "first" signals are subse- 
quently processed by a computer 16 as hereinafter described, to determine the angular position of a detected stator 
fault. 

[0046] The analogue signals from the notch filters 15 are also summed by a summing amplifier 17, the summed 
signal then being digitized. This summed signal must now be decomposed into its harmonic components. Two practical 

30 approaches to this task are possible. Given that the components of the axial leakage flux are mathematically predict- 
able, a set of filters, tuned to the frequencies of interest, can be used. The outputs of these filters can be monitored 
for changes which potentially indicate the onset of a fault condition. In practice, in view of the variety of motor designs 
which exist, construction of discrete filters tailored to the motor parameters is not a practical approach. However, by 
using digitally controlled analogue filters, it is possible to set up an array of filters which can be readily customized to 

35 the particular motor. For example, a single controllable filter may be sequentially stepped through the predetermined 
frequencies of interest. 

[0047] The preferred method, according to the present invention, of decomposing the signal into its harmonic com- 
ponents is to digitize the signal from the summing amplifier 17 and to transform the digital signal to the frequency 
domain by means of a fast Fourier transform (FFT). In this way, a further signal having a value corresponding to the 
40 magnitude of selected harmonic components of axial leakage flux is derived. The harmonic components are selected 
as being those which are significantly affected by the occurrence of a stator fault. 
[0048] The preferred method will now be described with reference to Figures 6 and 7. 

[0049] In order to detect the occurrence of a stator fault by monitoring changes in the axial leakage flux, it is first 
necessary to identify the frequencies of the harmonic components which are significantly affected by the occurrence. 

45 These can be calculated from the machine parameters and the conditions under which the motor is operating. 

[0050] Referring to Figure 6, the first thing to do is to input into the computer the relevant input parameters, specifically 
the number of poles of the motor and its nominal speed (BLOCK 21 ). From these parameters the harmonic components 
of axial leakage flux under faulted and no-fault conditions are computed (BLOCK 22). If the motor is on-load (BLOCK 
23), the rotor velocity has to be taken into account (BLOCK 24); in either case the results of the calculation are input 

50 into the initialization procedure (BLOCK 25). The summed outputs from the search coils are recorded at a first period 
t(n) and converted to the frequency domain, A(f,t(n)) (BLOCK 26) as previously described. The results are tested 
(BLOCK 27) and should confirm the calculations of (BLOCK 23). 

[0051] The next step in the procedure is to measure and record the summed output of induced emf on each of the 
search coils at a second period t(n+1 ), and converted to the frequency domain, A(f ,t(n+1 )) (BLOCK 29). L If a possible 
55 fault is detected by applying the test 



A(f,t(n+1)) » A(M(n)) 



(BLOCK 30) 
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then the summed output from the search coils is measured and recorded at a third period and converted to the frequency 
domain, -»A(f,t(n + 2)) (BLOCK 31). The possible fault is confirmed or ignored by testing whether 

5 A(f,t(n + 2)) » A(f ,t(n)) (BLOCK 32) 

[0052] If the result is positive, the fault is confirmed; otherwise it can be taken that the apparent fault was due to a 
temporary or spurious event. 

[0053] If the fault is confirmed, the rms value of induced emf in each of the search coils in the faulted condition of 
10 the motor is measured and recorded, (BLOCK 33). The emf value for the faulted condition is subtracted from the stored 
emf value for the healthy condition, for each search coil, and the modulus of the difference is taken (BLOCK 34). From 
the moduli value so obtained the angular position of the fault with respect to a datum position can be calculated, (BLOCK 
35), as will be described with reference to Figure 7. The occurrence of the stator fault, and its position in the stator 
winding, are then displayed (BLOCK 36). 
is [0054] Upon confirmation of the occurrence of a stator fault (BLOCK 32), the angular position of the fault is determined 
by using the algorithm illustrated in Figure 7. This algorithm involves the use of values derived from signals obtained 
after the occurrence of the fault, as will now be described, as well as the following values: 

R-, = radius of the circle on which the axes of the search coils lie, this value having been entered and 

20 stored initially; 

R 2 = mean radius of the end winding of the stator, this value having been entered and stored initially; 

EMF (healthy) = the values of said "first 0 signals derived from the search coils, these values having been stored 
25 prior to the occurrence of the fault. 

[0055] Referring to Figure 7, the number of search coils n was previously entered into the computer and stored 
(BLOCK 41). In the present example n = 4, and in the algorithm the individual search coils are identified by the letter 
x, where x = 1 , x = 2, x = 3 and x = 4, respectively. 

30 [0056] The signals from the search coils, previously referred to as the "first 0 signals, are examined in turn, (BLOCKS 
42, 43, 44). The "first" signal from each of the search coils corresponds to the axial leakage flux detected by the search 
coil, being derived from the voltage signal which is a measure of the rms value of the emf induced in the coil. The value 
of this signal, denoted by EMF (faulted, position x), is subtracted from the previously stored value of the signal EMF 
(healthy, position x), (BLOCK 45), and the modulus IB(x)l is taken (BLOCK 46). 

35 [0057] The computer now has in store all the necessary values and parameters from which the angular position of 
the stator fault can be computed, namely R-,, R 2 , B^ B 2 , B 3 , B 4 . From these values the following relationships are 
derived, (see BLOCK 47): 

40 ± ^ = cos' 1 [0.5(R 1 /R 2 +R 2 /R 1 )*(B 1 -B 2 )/(B 1 +B 2 )] 

± y 2 = cos" 1 [0.5(R 1 /R 2+ R 2 /R l ).(B 3 -B 4 )/(B 3+ B 4 )] 

45 since y = \|/ 2 + 90°, the required angular position of the fault \j/ 0 can be determined. In other words the quadrant in 
which the fault lies can be readily determined. 

[0058] Thus, considering the first pair of diametrically opposed search coils 1 , 2, one must determine whether EMF 
(fault, 1) is greater or less than EMF(fau!t, 2), (BLOCK 48). If it is greater, then the fault is at one of the positions ±\|/ 1 
from coil 1 , (BLOCK 49). If it is less, then the fault is at one of the positions ±\\f^ from coil 2, (BLOCK 50). In order to 
so resolve the ambiguity, one must now determine whether EMF(fault, 3) is greater or less than EMF(fault, 4), (BLOCK 
51). If it is greater, then the fault is at a position ±\|/ 2 from coil 3, (BLOCK 52), while if it is less the fault is at one of the 
positions ±\j/ 2 from coil 4, (BLOCK 53). The fault obviously lies within the quadrant defined by ±1^ and ±\|/ 2 , (BLOCK 
54); that is to say, the angular position \|/ 0 of the fault with respect to the datum position is given by 

55 

V 0 = V^ = V 2 +90°. 
This position is displayed along with indication that a fault has occurred. 
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Claims 

1. A method of detecting and locating stator faults in a dynamoelectric machine having a rotor and a stator with 
multiturn stator coils by monitoring changes in axial leakage flux during operation of the machine, comprising: 

identifying the frequencies of selected harmonic components of axial leakage flux that are to be significantly 
affected by the occurrence of a stator fault, 

detecting the axial leakage flux at each of at least four positions, including a datum position, distributed sym- 
metrically around the axis of the rotor adjacent to an end winding of the stator, 

deriving from the detected flux at each of said positions a first signal having a value corresponding to the 
magnitude of the detected flux, 

storing the values of said first signals derived prior to the occurrence of a stator fault, 

deriving from the detected axial leakage flux a further signal having a value corresponding to the magnitude 

of said selected harmonic components, 

monitoring said further signal and detecting a change therein thereby detecting the occurrence of a stator fault, 
storing the values of said first signals derived after the occurrence of the stator fault, and 
comparing the values of said first signals stored before and after the occurrence of a stator fault and storing 
the differences therebetween, 

characterized in that the method further comprises: 

determining from said comparison a value representing the angular position of the stator fault with respect to 
the datum position by: 

(i) storing a value representing the radial distance of said four positions from the axis of the stator, 

(ii) storing a value R 2 representing the mean radius of the stator end winding, 

(iii) computing values B-, and B 2 corresponding respectively to the moduli of said differences between the 
stored values of said first signals derived from the detected flux at one diametrically opposed pair of said 
flux detecting positions, 

(iv) computing values B 3 and B 4 corresponding respectively to the moduli of said differences between the 
stored values of said first signals derived from the detected flux at the other diametrically opposed pair of 
said flux detecting positions, 

(v) computing values ^ and \j/ 2 , where 

± Vi = cos* 1 [0.5(R 1 /R 2 +R 2 /R l )-(B 1 -B 2 )/(B 1 +B 2 )] 

± ¥2 = cos* 1 [0.5(R 1 /R 2+ R 2 /R 1 )*(B 3 -B 4 )/(B 3+ B 4 )], 

(vi) determining therefrom the angular position of the stator fault y 0 with respect to said datum position 
from the relationship 

¥ 0 = Vi =V 2 +90° 

and (vii) displaying the value \|/ 0 . 

2. A method according to claim 1 , wherein the storage of said first signals derived after the occurrence of a stator 
fault is initiated by the detection of the fault. 

3. A method according to claim 1 or 2, wherein the axial leakage flux is detected at each of four positions, including 
said datum position, distributed symmetrically around said axis. 

4. A method according to any one of claims 1 to 3, wherein the spectrum of axial leakage flux is analysed by 

deriving a respective voltage signal corresponding to the detected leakage flux at each of said plurality of 
positions, summing the derived voltage signals, digitising the sum of said voltage signals, and transforming the 
digitised signal to the frequency domain by means of a fast Fourier transform. 



EP 0 552 991 B1 ' 



5. Apparatus for detecting and locating stator faults in an induction motor (1 0) by monitoring changes in axial leakage 
flux during operation of the motor, comprising: 

means (1 3) for detecting the axial leakage flux during operation of the motor, said means including means for 
deriving a plurality of first signals each having a value corresponding to the detected flux at a respective one 
of at least four positions, including a datum position, distributed around the axis of the motor, 
means for deriving from the detected flux a further signal corresponding to the magnitude of the detected flux 
at selected harmonic frequencies, 

first data processing means (1 6) for monitoring said further signal and for detecting the occurrence of a stator 
fault by a change in the value thereof, and 

second data processing means (1 6) responsive to said first data processing means, said second data process- 
ing means including: 

(a) means for storing said first signals prior to the occurrence of a stator fault, and 

(b) means responsive to the detection of a stator fault for storing said first signals after the occurrence of 
the fault, 

characterized in that said second data processing means further includes: 

(c) means for computing from said stored values of said first signals the angular position of the stator fault 
with respect to said datum position, 

(d) said means for computing including: 

(d1 ) means for storing a value R., representing the radial distance of said four positions from the axis 
of the stator, 

(62) means for storing a value R 2 representing the mean radius of the stator end winding, 
(d3) means for computing values B A and B 2 corresponding respectively to the moduli of said differ- 
ences between the stored values of said first signals derived from the detected flux at one diametrically 
opposed pair of said flux detecting positions, 

(d4) means for computing values B 3 and B 4 corresponding respectively to the moduli of said differ- 
ences between the stored values of said first signals derived from the detected flux at the other dia- 
metrically opposed pair of said flux detecting positions, 
(d5) means for computing values and \|/ 2 , where 

± =cos* 1 [a5(R 1 /R 2 +R 2 /R l MB 1 -B 2 )/(B 1 +B 2 )] 

± V 2 = cos' 1 [0.5(R 1 /R 2+ R 2 /R 1 MB 3 -B 4 )/(B 3 +B 4 )], 

(d6) means for determining therefrom the angular position of the stator fault vj/ 0 with respect to said 
datum position from the relationship 

± ¥ 0 = ¥1 = V 2 + 90 ° 

and (d7) displaying the value \j/ 0 . 

6. Apparatus according to claim 5, wherein said means for detecting the axial leakage flux comprises four flux sensing 
coils (13) mounted symmetrically around the axis of the stator (11) adjacent to an end winding thereof. 



Patentanspruche 

1 . Verfahren zum Erfassen und Lokalisieren von Statorfehlern in einer dynamoelektrischen Maschine mit einem Rotor 
und einem Stator, mit Statorspulen mit mehreren Wicklungen, durch Uberwachen von Veranderungen des axialen 
Leckflusses wahrend des Betriebes der Maschine, mit den Merkmalen: 
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Identifizieren der Frequenzen ausgewahlter harmonischer Komponenten des axialen Leckflusses, die durch 
das Auftreten eines Statorfehlers in signifikanter Weise beeinfiuBt werden sollten, 

Erfassen des axialen Leckflusses an jeder von zumindest vier Positionen, einschlieBlich einer Ausgangspo- 
sition, welche symmetrisch urn die Achse des Rotors in der Nahe einer Entwicklung des Stators verteilt sind, 

Ableiten eines ersten Signals aus dem erfaBten FluB an jeder der Positionen, wobei das Signal einen Wert 
hat, welcher der GroBe des erfaBten Flusses entspricht, 

Speichern der Werte der ersten Signale, die vor dem Auftreten eines Statorfehlers abgeleitet wurden, 

Ableiten eines weiteren Signals aus dem erfaBten axialen LeckfluB, wobei das Signal einen Wert hat, welcher 
der GroBe der ausgewahlten harmonischen Komponenten entspricht, 

Qberwachen des weiteren Signals und Erfassen einer Veranderung desselben, urn dadurch das Auftreten 
eines Statorfehlers zu erfassen, 

Speichern der Werte der nach dem Auftreten des Statorfehlers abgeleiteten ersten Signale, und 

Vergleichen der Werte der ersten Signale, die vor und nach dem Auftreten eines Statorfehlers gespeichert 
wurden und Speichern der Unterschiede dazwischen, 

dadurch gekennzeichnet, daB das Verfahren weiterhin aufweist: 

aus dem Vergleich Bestimmen eines Wertes, welcher die Winkelposition eines Statorfehlers bezuglich der 
Ausgangsposition wiedergibt, durch: 

(i) Speichern eines Wertes welcher den radialen Abstand der vier Positionen von der Achse des Stators 
wiedergibt, 

(ii) Speichern eines Wertes R 2 , welcher den mittleren Radius der Statorentwicklung wiedergibt, 

(iii) Berechnen von Werten B n und B 2 , die jeweils den Amplituden der Unterschiede zwischen den gespei- 
cherten Werten der ersten Signale entsprechen, welche von dem erfaBten FluB an einem diametral ge- 
genuberliegenden Paar der FluBerfassungspositionen abgeleitet wurden, 

(iv) Berechnen von Werten B 3 und B4, welche jeweils den Amplituden der Unterschiede zwischen den 
gespeicherten Werten der ersten Signale entsprechen, die von dem erfaBten FluB an dem andern dia- 
metral gegenuberliegenden Paar von FluBerfassungspositionen abgeleitet wurden, 

(v) Berechnen von Werten ^ und \j/ 2 , wobei 

+ = cos" 1 [0,5(R l /R 2 +R 2 /R 1 HB 1 -B 2 )/(B 1 +B 2 )] 

± ¥2 = cos- 1 [0,5(R 1 /R 2+ R 2 /R 1 ).(B 3 -B 4 )/(B 3+ B 4 )] 

(vi) daraus Bestimmen der Winkelposition des Statorfehlers \|/ 0 bezuglich der Ausgangsposition gemaB 
der Beziehung 

V 0 = Vi = V 2 + 90<> 

und (vii) Anzeigen des Wertes y 0 . 

Verfahren nach Anspruch 1 , wobei das Speichern der ersten Signale, welche nach dem Auftreten eines Stator- 
fehlers abgeleitet wurden, durch die Erfassung des Fehlers ausgelost wird. 
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Verfahren nach Anspruch 1 oder 2, wobei der axiale LeckfluB an jeder der vier Posltionen erfaBt wird, einschlieBlich 
der Ausgangsposition, welche symmetrisch urn die Achse verteilt sind. 

Verfahren nach einem der Anspruche 1 bis 3, wobei das Spektrum des axialen Leckflusses analysiert wird durch 
Ableiten jeweils eines Spannungssignals, welches dem erfaBten LeckfluB an jeder der Mehrzahl von Posi- 
tionen entspricht, Aufsummieren der abgeleiteten Spannungssignale, Digitalisieren der Summe der Spannungs- 
signale und Transformieren des digitalisierten Signals in den Frequenzraum bzw. die Frequenzdomane mit Hilfe 
einer schnellen Fourier-Transformation. 

Vorrichtung zum Erfassen und Lokalisieren von Statorfehlern in einem Induktionsmotor (10) durch Uberwachen 
von Veranderungen in dem axialen LeckfluB wahrend des Betriebes des Motors, mit: 

einer Einrichtung (1 3) zum Erfassen des axialen Leckflusses wahrend des Motorbetriebes, wobei die Einrich- 
tung Einrichtungen fur das Ableiten einer Mehrzahl erster Signale aufweist, die jeweils einen Wert haben, 
welcher dem erfaBten FluB an jeweils einer von zumindest vier Positionen entspricht, und zwar einschlieBlich 
einer Ausgangsposition, welche urn die Achse des Motors verteilt sind, 

Einrichtungen fur das Ableiten eines weiteren Signals aus dem erfaBten FluB, welches der GroBe des erfaBten 
Flusses bei ausgewahlten harmonischen Frequenzen entspricht, 

einer ersten Datenverarbeitungseinrichtung (16) zum Uberwachen des weiteren Signals und zum Erfassen 
des Auftretens eines Statorfehlers durch eine Veranderung in dessen Wert, und 

einer zweiten Datenverarbeitungseinrichtung (16), die auf die erste Datenverarbeitungseinrichtung reagiert, 
wobei die zweite Datenverarbeitungseinrichtung aufweist: 

(a) Einrichtungen zum Speichern der ersten Signale vor dem Auftreten eines Statorfehlers, und 

(b) Einrichtungen, welche auf die Erfassung eines Statorfehlers ansprechen, urn die ersten Signale nach 
dem Auftreten des Fehlers zu speichern, 

dadurch gekennzeichnet, daB die zweite Datenverarbeitungseinrichtung weiterhin aufweist: 

(c) Einrichtungen zum Berechnen der Winkelposition des Statorfehlers bezuglich der Ausgangsposition 
aus den gespeicherten Werten der ersten Signale, 

(d) wobei die Einrichtung fur das Berechnen aufweist: 

(d1) Einrichtungen fur das Speichern eines Wertes R v welcher den radialen Abstand der vier Posi- 
tionen von der Achse des Stators wiedergibt, 

(62) eine Einrichtung fur das Speichern eines Wertes R 2 , welcher den mittleren Radius der Stator- 
entwicklung wiedergibt, 

(d3) eine Einrichtung fur das Berechnen von Werten und B 2 , welche jeweils den Amplituden der 
Unterschiede zwischen den gespeicherten Werten der ersten Signale entsprechen, die aus dem er- 
faBten FluB an einem diametral gegenuberliegenden Paar der FluBerfassungspositionen erfaBt wur- 
den, 

(d4) eine Einrichtung fur das Berechnen von Werten B 3 und B 4 , welche jeweils den Amplituden der 
Unterschiede zwischen den gespeicherten Werten der ersten Signale entsprechen, die aus dem er- 
faBten FluB an dem anderen diametral gegenuberliegenden Paar von FluBerfassungspositionen ab- 
geleitet wurden, 

(d5) eine Einrichtung fur das Berechnen von Werten vj^ und ij/ 2 , wobei 



± =cos" 1 [0,5(R 1 /R 2 +R 2 /R 1 MB 1 -B 2 )/(B 1 +B 2 )] 
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± y 2 = cos°[0,5(R 1 /R 2+ R 2 /R 1 WB3-B 4 )/(B 3+ B 4 )] l 

(d6) eine Einrichtung, um daraus die Winkelposition des Statorfehlers \j/ 0 bezuglich der Ausgangspo- 
sition zu bestimmen, und zwar gemaB der Beziehung 

± V 0 = Vi = ¥ 2 + 90 ° 

und (d7) Anzeigen des Wertes y 0 . 

6. Vorrichtung nach Anspruch 5, wobei die Einrichtung fur das Erfassen des axialen Leckflusses vier FluBsensor- 
spulen (13) aufweist, die symmetrisch um die Achse des Stators (11) in der Nahe einer Endwicklung desselben 
montiert sind. 



Revendications 

1 . Proc6de de detection et de localisation de defaillances du stator dans une machine dynamo6lectrique, comportant 
un rotor et un stator k bobines statoriques multispires, par une surveillance des modifications du flux axial de fuite 
pendant le fonctionnement de la machine, comprenant : 

Identification des frequences de composantes harmoniques s6lectionn6es du flux axial de fuite, qui sont 
affectees de maniere significative par Papparition d'une defaillance du stator, 

la detection du flux axial de fuite au niveau de chacune d'au moins quatre positions, y compris une position 
de reference, reparties symetriquement autour de Paxe du rotor, de maoi&re adjacente a un enroulement 
d'extremite du stator, 

I'obtention, k partir du flux detecte au niveau de chacune desdites positions, d'un premier signal ayant une 
valeur qui correspond k la grandeur du flux detects, 

la memorisation des valeurs desdits premiers signaux obtenus avant Papparition d'une defaillance du stator, 
I'obtention, k partir du flux axial de fuite detects, d'un autre signal ayant une valeur qui correspond k Pamplitude 
desdites composantes harmoniques s^lectionnees, 

la surveillance dudit autre signal et la detection d'une modification dans celui-ci, pour ainsi detecter Papparition 
d'une defaillance du stator, 

la memorisation des valeurs desdits premiers signaux, obtenus aprfcs Papparition de la defaillance du stator, et 
la comparaison des valeurs desdits premiers signaux, memorises avant et apres Papparition d'une defaillance 
du stator, et la memorisation des differences entre ceux-ci, 

caracterise en ce que le proc6de comprend en outre : 

la determination, k partir de ladite comparaison, d'une valeur repr6sentant la position angulaire de la defaillan- 
ce du stator par rapport k la position de reference : 

(i) en enregistrant une valeur R v representant la distance radiale desdites quatre positions par rapport k 
Paxe du stator, 

(ii) en enregistrant une valeur R^ representant le rayon moyen de Penroulement d'extremite du stator, 

(iii) en calculant des valeurs et B 2 , correspondant respectivement aux modules desdites differences 
entre les valeurs memorisees desdits premiers signaux, obtenus k partir du flux detecte au niveau d'une 
paire desdites positions de detection de flux, diametralement opposees, 

(iv) en calculant des valeurs B3 et B 4 , correspondant respectivement aux modules desdites differences 
entre les valeurs m6moris6es desdits premiers signaux, obtenus k partir du flux detecte au niveau de 
{'autre paire desdites positions de detection de flux, diametralement oppos6es, 

(v) en calculant des valeurs et \y 2 , ou 



± Vl = cos" 1 [O^R/R,, + R 2 /R 1 ) . (B 1 - B 2 )/(B 1 + B 2 )] 



EP 0 552 991 B1 



± y 2 = cos -1 [0,5(FyR 2 + R 2 /R 1 ) . (B 3 - B 4 )/(B 3 + B 4 )], 

(vi) en determinant, sur cette base, la position angulaire de la defaillance du stator y 0 par rapport k ladite 
position de reference, k partir de la relation 

¥0 = ^1 = V 2 + 90 °» 

et 

(vii) en affichant la valeur y 0 . 

Proced6 selon la revendication 1 , dans lequel la memorisation desdits premiers signaux obtenus apres I'apparition 
d'une defaillance du stator est Ianc6e par la detection de la defaillance. 

Proc6de selon la revendication 1 ou 2, dans lequel le flux axial de fuite est detecte au niveau de chacune de quatre 
positions, y compris ladite position de reference, ^parties symetriquement autour dudit axe. 

Procede selon Tune quelconque des revendications 1 k 3, dans lequel le spectre du flux axial de fuite est analyse 
par 

i'obtention d'un signal de tension respectif, correspondant au flux de fuite detecte au niveau de chacune de 
ladite plurality de positions, I'addition des signaux de tension obtenus, la num6risation de la somme desdits signaux 
de tension, et la transformation du signal numeris6 en domaine frequenciel au moyen d'une transform^ de Fourier 
rapide. 

Appareil de detection et de localisation de defaillances du stator dans un moteur k induction (10) par une sur- 
veillance des modifications du flux axial de fuite pendant le fonctionnement du moteur, comprenant : 

un moyen (13) servant k detecter le flux axial de fuite pendant le fonctionnement du moteur, ledit moyen 
comportant un moyen permettant d'obtenir une pluralite de premiers signaux, ayant chacun une valeur cor- 
respondant au flux detects au niveau de Tune d'au moins quatre positions, y compris une position de reference, 
reparties autour de I'axe du moteur, 

un moyen permettant d'obtenir, k partir du flux detecte, un autre signal correspondant a la grandeur du flux 
d6tect6 k des composantes harmoniques seiectionnees, 

un premier moyen de traitement de donnees (16) prevu pour surveiller ledit autre signal et pour detecter 
I'apparition d'une defaillance du stator, par une modification de la valeur de ce signal, et 
un second moyen de traitement de donnees (1 6), rSagissant audit premier de moyen de traitement de donn6es, 
ledit second moyen de traitement de donnees comprenant : 

(a) un moyen pour memoriser lesdits premiers signaux, avant I'apparition d'une defaillance du stator, et 

(b) un moyen reagissant k la detection d'une defaillance du stator pour memoriser lesdits premiers signaux 
apres I'apparition de la defaillance, 

caracteris6 en ce que ledit second moyen de traitement de donnees comprend en outre : 

(c) un moyen pour calculer, k partir desdites valeurs memorisees desdits premiers signaux, la position 
angulaire de la defaillance du stator par rapport k ladite position de reference, 

(d) ledit moyen pr6vu pour le calcul comprenant : 

(d1 ) un moyen pour enregistrer une valeur representant la distance radiale desdites quatre posi- 
tions par rapport k I'axe du stator, 

(d2) un moyen pour enregistrer une valeur R 2 , representant le rayon moyen de Penroulement d'ex- 
tremite du stator, 

(d3) un moyen pour calculer des valeurs et B 2 , correspondant respectivement aux modules des- 
dites differences entre les valeurs memoris6es desdits premiers signaux, obtenus k partir du flux 
detecte au niveau d'une paire desdites positions de detection de flux, diam6tralement opposees, 
(d4) un moyen pour calculer des valeurs B 3 et B 4 , correspondant respectivement aux modules des- 
dites differences entre les valeurs memorisees desdits premiers signaux, obtenus k partir du flux 
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detects au niveau de Pautre paire desdites positions de detection de flux, diam6tralement oppos6es, 
(d5) un moyen pour calculer des valeurs ^ et \j/ 2 , ou 

± = cos' 1 [0,5(R 1 /R 2 + Rg/F^) . (B 1 - B 2 )/(B 1 + B 2 )] 

± V 2 = cos -1 [O.SJR/Rg + R 2 /R t ) . (B 3 - B 4 )/(B 3 + B 4 )], 

(d6) un moyen permettant de determiner, sur cette base, la position angulaire de la defaillance du 
stator y 0 P ar rapport & ladite position de reference, k partir de la relation 

± Vo = \|/ 1 = ¥2 + 90°, 

et 

(d7) Paffichage de la valeur y 0 . 

Appareil selon la revendication 5, dans lequel ledit moyen servant & detecter le flux axial de fuite comprend quatre 
bobines detectrices de flux (13) montees de manure sym&rique autour de Paxe du stator (11), de maniere adja- 
cente & un enroulement d'extr6mit§ de celui-ci. 
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((Bi-B2)/(Bi*B2))l 
1^2 = COS^ I0.5(R 1 /R 2 *R 2 /Ri) 
((B 3 -B4)/(B3*B 4 ))] 
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